Gravitating fluids with Lie symmetries 



A M Msomif t, K S Govinderfand S D Maharajf 

f Astrophysics and Cosmology Research Unit, School of Mathematical Sciences, 
University of KwaZulu-Natal, Private Bag X54001, Durban 4000, South Africa 
Email: maharaj@ukzn.ac.za 

Abstract. We analyse the underlying nonlinear partial differential equation which 
arises in the study of gravitating flat fluid plates of embedding class one. Our interest 
in this equation lies in discussing new solutions that can be found by means of Lie 
point symmetries. The method utilised reduces the partial differential equation to an 
ordinary differential equation according to the Lie symmetry admitted. We show that 
a class of solutions found previously can be characterised by a particular Lie generator. 
Several new families of solutions are found explicitly. In particular we find the relevant 
ordinary differential equation for all one-dimensional optimal subgroups; in several 
cases the ordinary differential equation can be solved in general. We are in a position 
to characterise particular solutions with a linear barotropic equation of state. 
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1. Introduction 

The local isometric embedding of four-dimensional Riemannian manifolds M4 in 
higher dimensional flat pseudo-Euchdean spaces En{N < 10) is important for several 
applications in general relativity. For the basic theory and general results pertinent 
to embeddings the reader is referred to Stephani et al (2003). The invariance of the 
embedding class naturally generates a classification scheme for all solutions of the field 
equations in terms of their embedding class. The embedding class p is the minimum 
number of extra dimensions of the Riemannian manifold M4, ie. p — N — 4. Exact 
solutions have been found by the method of embedding in particular spacetimes for 
simple cases of low embedding class. Some of these exact solutions may not be easily 
found using other methods and techniques. For example the embedding method has been 
utihsed to find all conformally fiat perfect fiuid solutions, in embedding class p = 1, of 
Einstein's field equations (Krasinki 1997, Stephani 1967a, Stephani 1967b). We point 
out that embedding of four-dimensional Riemannian manifolds in higher dimensional 
spacetimes with arbitrary Ricci tensors has been investigated by several authors. The 
physical motivation here is to understand the nature of physics in higher dimensions; 
the modern view is that the Riemannian manifold M4 is a hypersurface in the higher 
dimensional bulk in the brane world scenario and other higher dimensional themes 
(Dahia and Romero 2002a, Dahia and Romero 2002b, Dahia et al 2008). 

Gupta and Sharma (1996) have generated a relativistic model in higher dimensions 
describing gravitating fiuid plates. Advantages of this model are that it is easy to 
interpret the physical features using embedding in higher dimensions and the underlying 
differential equation governing the gravitational dynamics is tractable. This model is 
expanding and not conformally fiat. A plane symmetric metric in four-dimensional 
spacetimes M4 given by 

ds^ = -dr' - t\d9^ + e^d(t>^) + (1 + 2V)dt^ + 2V'drdt (1) 

is embedded in the five-dimensional pseudo- Euclidean space £"5 with metric 

ds^ = -{dz^ - {dz^Y - {dz^ + {dz^ - (d^5)2. (2) 

This embedding is achieved by setting 

z^ — to cos (j) (3a) 
z'^^tesincj) (36) 

(3c) 
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+ V {3d) 

(3e) 

where V = V{r, t) is an arbitrary function. 

Consequently this model has embedding class p = 1 which allows both conformally 
flat and nonconformally flat fluid distributions. The solutions admitted may be geodesic 
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or accelerating. For a nonzero conformal (Weyl) tensor a partial differential equation 
has to be satisfied. This pivotal equation governs the evolution of the system and a 
particular class of solutions was identified by Gupta and Sharma (1996) by inspection. 
A detailed analysis of the pivotal equation shows that other classes of solution are 
possible which contain the Gupta and Sharma (1996) models as a special case. 

Our intention is to systematically study the pivotal equation and to obtain a deeper 
insight into the nature of solutions permitted using the Lie analysis of differential 
equations. In §2 we discuss the fundamental partial differential equation that governs 
the gravitational behaviour of the model, and present known solutions. An outline of 
the basic features of the Lie symmetry analysis is given in §3. We regain the Gupta and 
Sharma (1996) models using the relevant Lie generator in §4. In §5 we systematically 
study the group invariant solutions admitted by the fundamental equation. The partial 
differential equation is reduced to an ordinary differential equation for each element of 
the optimal system. The integrability of the ordinary differential equation is considered 
in each case and exact solutions are identified. The physical features of the model are 
discussed in §6 and equations of state are found for particular solutions. Some brief 
concluding comments are made in §7. 



2. The model 



The embedding of the four-dimensional Riemannian metric (1) into the five-dimensional 
fiat metric (2) leads to a differential equation that is central to the model. Gupta and 
Sharma (1996) show that the pivotal equation is 



-VttVrr + Vrt' + ^ [(1 + 2K)K-r- - Vu - 2KtK] 



+ 



t2 



l + 2Vt + Vr 







(4) 



where subscripts denotes partial differentiation. This is a nonlinear equation in V and 
difficult to solve. We need to explicitly solve (4) to describe the gravitational dynamics. 

To demonstrate a class of solutions to (4), Gupta and Sharma (1996) made the 
following assumption 

7(r) 



+ h{t) + Ci 



where C and Ci are arbitrary constants. Then (4) reduces to the separable form 



4/- 



ACfr 



- - thu + 2ht 



a 



(5) 



(6) 



C/rr + 1 

where a is the constant of separability. It is possible to solve the equation (6) in terms 
of t explicitly as 

a t ait^ 
Tt~2C 3~ 



+ 0-2 



(7) 
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and to provide four solutions to the equation in terms of r, viz 

1 la 

^-^^'^'' + ^■ + 1 

1 rv 

/=^(r + /3f+4 (8c) 

/ = coshX + - (8d) 

where we have set X — \/2mr + mo and cui, Q!2, /5, m and mo are arbitrary constants. 

Equations (7) and (8a)-(8cO are then combined to provide solutions to the original 
equation (4), viz 

1 It kt^ 

'' = 2^t''^'l^2;;^t--2-T^'^ 

V=Ytir + (3r-'--^ + h (9c) 

V^^coshX-^--- — + ki (9d) 
2mH 2mH 2 3 ^ ' 

where k = Cai and ki = Ca2- Thus the assumption (5) leads to a simple class of 

solutions (9a)-(9(i) which are written in terms of elementary functions. As an aside 

we observe that a does not appear in the solutions {9a)-{9d). Thus the constant of 

separability a can be taken to be zero, and C can be taken to be unity with no loss of 

generality. Indeed as we shall demonstrate later, our Lie analysis obviates the need for 

the introduction of a and C. 

We will show that, while (5) is an ad hoc assumption, the reason for its feasibility 

lies in the group theoretic properties of (4). Utilising the full group properties of (4) we 

can provide further solutions to complement (9a)-(9(i). 

3. Lie analysis 

The basic feature of Lie analysis requires the determination of the one-parameter (e) 
Lie group of transformations 

t ^t + sr(t, r, u) + 0{e^) 

r ^r + e^{t,r,u) + 0{e'^) (10) 
U — U + er]{t, r, u) + 0{e^) 

that leaves the solution set of a differential equation invariant. (The full details can be 
found in a number of excellent texts (Bluman and Kumei 1989, Olver 1993).) In order 
to obtain (10) we need to determine their "generator" 
d d d 



Gravitating fluids with Lie symmetries 



5 



(also called a symmetry of the differential equation) which is a set of vector fields. 

The determination of these generators is a straight forward, albeit tedious process. 
Fortunately, a number of computer algebra packages are available to aid the practitioner 
(Hereman 1994). While some modern packages have been developed (Dimas and 
Tsoubelis 2005, Cheviakov 2007), we have found the package PROGRAM LIE (Head 1993) 
to be the most useful in practice. Indeed it is quite remarkable how accomphshed such 
an old package is - it often outperforms its modern counterparts! 

Utihsing PROGRAM LIE, we can demonstrate that (4) admits the following Lie point 
symmetries / vector fields : 

G, = I, (12a) 

G, = ^ (126) 

Gs = f^ (12c) 

«' = '^ + '-|7 (12<i) 

G« = r|: + (t + 2V)^ (12i) 
with the nonzero Lie bracket relationships 

[Gi, Gs] = —Gi [Gi, Gq] = 2Gi 

[G2, G4] = Gi [G2, Gq] = G2 

[G3, G5] = — 4G3 [G3, Gq] = 2G3 

[(74, Gs] = — G4 [G4, Gq] = G4 

for the given fields. As a result, the symmetries (12a)-(12/) form a six-dimensional 
indecomposable solvable Lie algebra, L (Rand et al, 1988). While L is not nilpotent, its 
first derived Lie subalgebra L^^^ —< Gi, G2, G3, G4 > is nilpotent and also represents the 
nilradical of L. Further information about such Lie algebras can be found in (Turkowski, 
1990). 

4. Known solutions 

It is possible to demonstrate that the solutions (9a)-(9d) are a natural consequence of 
a subset of these symmetries. If we take the combination 

G = ciGi + C2G2 + G5 

= 4 + (c.-( + crf'-V)|7 (14) 

the partial differential equation (4) is reduced to the ordinary differential equation 

UUrr -Ur^ + U = (15) 

as the essential equation governing the gravitational dynamics where V = c^ — ^ + ^^ + j. 



(13) 
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On comparing (6) and (15) we can identify the function / with U. We are now 
in a position to make a number of comments relating to the underlying assumption (5) 
in the Gupta and Sharma (1996) solutions. We observe that the t-dependence arises 
naturally because of the choice of the symmetry G. Thus the temporal dependence is 
not arbitrary as suggested by the function h{t) in the choice (5). It is not necessary to 
solve any differential equation to obtain the form of h{t) given by (7). In addition, the 
solutions of (15) are the same as {8a)-{8d) with a — and C — 1. As stated earlier, 
the final solutions (9a)-(9d) do not include these constants. Thus the Lie symmetry G 
leads directly to the canonical form of the solution to (4) without the need to introduce 
spurious arbitrary functions and parameters. 

5. Group inveiriant solutions 

We now seek to utilise the Lie point symmetries in a systematic manner to generate 
new solutions. These new solutions are termed group invariant solutions as they will be 
invariant under the group generated by the symmetry used to find them. The advantage 
of using Lie point symmetries is that we are guaranteed that the variable combinations 
obtained will always result in an equation in the new variables - no further "consistency" 
conditions are needed. As we arc dealing with a 1 + 1 partial differential equation here, 
wc will always be able to find an ordinary differential equation in the new variables 
defined by the symmetries. 

5.1. The optimal system 

Given that equation (4) has the six symmetries (12a)-(12/), we can find group invariant 
solutions using each symmetry individually, or any linear combination of symmetries. 
However, taking all possible combinations into account is overly excessive. It turns out 
(Olver, 1993), that one only need consider a subspace of this vector space. We use the 
subalgebraic structure of the symmetries (12a)-(12/) of the system (4) to construct an 
optimal system of one-dimensional subgroups. Such an optimal system of subgroups 
is determined by classifying the orbits of the infinitesimal adjoint representation of a 
Lie group on its Lie algebra obtained by using its infinitesimal generators. All group 
invariant solutions can be transformed to those obtained via this optimal system. The 
process is algorithmic and can be found in (Olver 1993). Here we only summarise the 
final results. 

In order to obtain group invariant solutions of (4) explicitly, the optimal system 
yields only the following symmetry combinations 




d 



(16 a) 



(166) 
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Table 1. One symmetry: generators and ODES. 



Generator 


Invariants 


ODE 


G2 


y = t 


yUyy -2Uy-l=Q 




V = U{t) 




G3 


y = t 


No ODE exists 




U{t) = r 




G4 


y = t 


yUyy -2Uy-l=Q 




y = i + u{t) 




G5 


y = r 


UUyy - {Uyf + [/ = 




V = -l + lU{r) 




Table 2. One symmetry: PDE solutions. 


Generator 


Solution to PDE 




G2 


V = A+ \Bt^ - 1 






No solution to pde 





Gi V ='^+ A+\Bt^ - 

Gs l/ = -| + i [2^(sin(x/2mr + mo) + l)] 

v = -l 

V = -l + J [2^(cosli(x/2mr + mo) - 1)] 



^'''^ 

G, + G4 = (l + t)|: + rA (16/) 

^^ + ^^ = 1 + 4-^^ + ^)^ » 

G, + G, = {t' + r)^ + t^ (16/.) 

aCs + Ge = r— + (a + + [(1 - a)V - at] — (16i) 
G2-Gs + G, = {t + l)— + {r- 1')— (16j) 

G2 + G3 + G^^it+1)— + (t^ + r) — . (16A;) 

All solutions of (4) which are obtained via other combinations of point symmetries can 
be transformed into the solutions obtained from the combinations above. Here, we have 
also taken into account the fact that (4) is invariant under the following involutions: 
> — t,r— > — r and V ^ —V and so were able to restrict the optimal system further. 
It is clear that the optimal system consists of single elements of the Lie algebra, 
combinations of two elements and combinations of three elements only. We divide our 
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Table 3. Two symmetries: generators and ODES. 


Generator 


Invariants 


ODE 




t = v 


wf7„„ - 2[/„ - 4w'5 - 1 = 

y^yy ^y y 




V = rt^ + U{t) 




G2 + G4 


t = y 


{l + 2y){-l-2Uy + yUyy)=0 








G3 + G4 


t = y 


2 + y^+AUy- 2yUyy = 




V = rt'' + ^ + U{t) 




G2+G5 


r = \nt + y 


Uyy + 4C/ + 4UUyy + 5Uy 




V = -^ + ^U{r-lnt) 


-4C/2 + UyUyy = 


0G5 + Ge 


r = 


-il + 2a)^U^ + y^Uyy 




V = -l+ t(i-«)/(i+«)f/(ri-^/(^+")) 


+2{-l + a){l + 2a)U{l + Uyy) 






+yUy{2 + 5a + aUyy) = 



Table 4. Two symmetries: VDE solutions. 

Generator Solution to pde 

G2 + G3 V = rt^+A+ ^Bt^ + y - I 

G3 + G4 V = rt^ + ^ + A+ iBt^ + M - I 
G2 + G5 No solution to pde 

aG, + G, y = -1 + + 1 (e^ - if )^/^ + 4^)) 

^" = -5) 

_ _ t _ , C^r" , VC^r^-4t^ _ Cr^VC^r*-4t^ , 
^ ~ 2 2t 12*3" T 3(^j J2t3 I" 

(a = l) 



discussion of tiie solutions based on tiiis separation. 
5.2. One generator 

In this section we generate solutions to the master equation (4) when it admits the 
individual Lie symmetries (^2,^3,(^4 and G5. We do not consider the generators Gi 
and Gq as they do not appear in the optimal system (16a)-(16A;). The procedure of 
generating the invariants, the resultant ordinary differential and finally the solution to 
the partial differential equation is standard. Consequently we do not provide the details 
of the calculations, instead the relevant results are collated in tabular form. 

In table 1 we present the ODES that are generated when a single Lie symmetry 
generator is present. In table 2 we give solutions to the partial differential equation (4) 
when it is invariant under a single generator from the optimal system (16a)-(16/^). For 
the generator G^ we do not obtain an invariant involving V and so it is not possible to 
generate an ordinary differential equation. Thus there is no solution possible invariant 
under Gs alone. For the generators G2 and G4 it is possible to solve the resulting 
ordinary differential equations and obtain explicit forms for the function V given in table 
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Generator Invariants ode 



G2 + G3+G4 t = y 2/(5 + y(9 + 2/(7 + y(2 + 2/2(2 + j;)2)))) 

^ = 2(f^ + (ffe) + U{t) -(1 + y?{l + 22/)(-2C/, + yUyy) = -1 

G2-G3 + G4 t = y y(5 + 2/(9 + y(7 + 2/(2 + y2(2 + y)2)))) 

V=^-^ + Uit) -(1 + yf{l + 2y){-2Uy + yUyy) = -1 



Table 6. Three symmetries: PDE solutions. 
Generator Solution to pde 



62 + G3 + G4 V - 2^Y+t) + XTTfJ + ^+ |(f~^ + T + ^ + TTt) 
+if'((l + S|)-|(l + 8t3)log(l + 2i)) ^ 

G2 - G3 + G4 V = 2(l+4y + ^fq^ + ^+ |(f-^ + T + ^ + TTt) 
+it3((i + B|)-|(l + 8t3)log(l + 2i)) 



2. These are new solutions to equation (4) which have not been obtained previously. 
For the generator G5 we have obtained the ordinary differential equation 

UUyy -Uy^ + U = ( 1 7) 

which is of the same form as (15) in §4. However it is important to observe that the 
characteristics here are different from §4. Consequently the solutions V generated by 
the Lie symmetry G5, and listed in table 2, comprise a new class of exact solutions to 
equation (4). 

5.3. Two generators 

We now consider the combinations of two generators which arise in the optimal system 
(16a)-(16A;). Table 3 contains the invariants and reduced ode. Table 4 lists the analytic 
solutions of the pde, containing all cases that we were able to obtain exphcit general 
solutions. 

For the generator G2 + G5, we have not been able to find a solution to the resultant 
ODE which is highly nonlinear. For the generator aG^+Ge the resultant ODE is nonlinear 
and difficult to solve in general. However particular solutions can be found for the special 
parameter values a — and a — 1. It is unlikely that the ODE will yield closed form 
solutions for other values of a. 

Thus we have generated new analytic solutions to (4) via combinations of two 
symmetries of (4) except for one case corresponding to the generator G2 + G5. 

5.4. Three generators 



It finally remains to consider the two combinations of three generators G2 + G3 + G4 
and G2 — Gs + G4which appear in the optimal system (16a)-(16A;). In both cases it is 
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possible to generate the invariants and reduced ODE which are presented in table 5. The 
explicit sohitions of the PDE are given in table 6. Both functions obtained in table 6 are 
new solutions to equation (4). 



6. Physical features 

We briefly describe the behaviour of the thermodynamic variables corresponding to the 
spacetime (1). The energy density and the pressure are given by 

= ^ + ^i^ttVrr - Vrt') (18a) 

Snp = ^ (186) 
respectively, where we have set 

P=l + 2Vt + Vr^. (19) 

For many applications in cosmology it is necessary that there exist barotropic equations 
of state in the form p = p{p) (Stephani et al 2003). We find that for the case of a single 
generator of the optimal system (16a)-(16A;) considered in §5.2 there exists a linear 
equation of state. The relevant equations of state are presented in Table 7. 

The equations of state p = p and p = \p were identifled by Gupta and Sharma 
(1996) for their class of solutions. We have demonstrated that their result follows 
because of the existence of the symmetry G5. The Lie symmetry G2 produces a new 
solution with equation of state p = p- The generator G4 gives another new solution with 
the linear equation of state p — y+Ij^P- Such linear equations of state are of importance 
in relativistic stellar structures and arise in models of quark stars (Komathiraj and 
Maharaj 2007, Mak and Harko 2004, Sharma and Maharaj 2007, Witten 1984). Also, 
in the modelling of anisotropic relativistic matter in the presence of the electromagnetic 
fleld for strange stars and matter distributions, we need a linear barotropic equation of 
state (Lobo 2006, Thirukkanesh and Maharaj 2008). 

For the generators considered in §5.3 and §5.4 there are no simple barotropic 
equations of state connecting the energy density and the pressure. However it is possible 
to describe the thermodynamical behaviour graphically. As an example we consider the 
solution corresponding to the generator G2 + G^. The energy density is 

_ 47r(5 + 3r) + 3(-3 + 2%^ 
^ ~ 167tH\2B + 6r + 3^4)2 

and pressure is 

^~ 8ivt\2B + 6r + 3t^)' 

Clearly there is no barotropic equation of state in this case. The behaviour of energy 
density has been plotted in figure 1 and the pressure is represented in figure 2. We have 
generated these plots with the help of Mathematica (Wolfram 1999). It is clear that there 
exist regions of spacetime in which p and p are well behaved, remaining finite, continuous 
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Generator Solution to pde Equation of state 



G2 


V = A+ \Bt^ - 1 


p 


G3 


No solution to pde 


No equation of state 








G5 




P = P 






P=\P 



and bounded. It is then viable to study the behaviour of the thermodynamical quantities 
such as the temperature over this region. We point out that plots for p and p for the 
other combinations of generators in the optimal system have similar behaviour. 




Figure 1. Graph of energy density for G2 + Ga. 




Figure 2. Graph of pressure for G2 + G3. 
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7. Conclusion 

A variety of new exact solutions of the governing equation (4), using the Lie method of 
infinitesimal generators, have been obtained. Previously known solutions where shown 
to be characterized by a particular Lie generator and are contained, as a special case, 
in our new family of solutions. We considered each element in the optimal system of 
one-dimensional subgroups and reduced the master equation to an ordinary differential 
equation. We were in a position to solve the resulting equations and obtain several new 
solutions for the gravitating model. A pleasing feature of our analysis is that several 
models generated admit a linear barotropic equation of state. The pivotal equation has 
six Lie point symmetries with eight nonzero Lie bracket relationships which generates 
the optimal system. It is this geometric structure which has enabled us to show that 
equation (4) has a rich structure. 
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